Cosmic ray anomalies and DAMA experiment in an Extended Seesaw Model 
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We show that an extended seesaw model proposed to achieve low scale leptogenesis can resolve the 
excess positron and electron fluxes observed from PAMELA, ATIC and/or Fermi-LAT, and simulta- 
neously accommodate the recent DAMA experimental results. In this approach, an extra vector-like 
singlet neutrino S and a singlet scalar field ip, which are doubly coexisting dark matter candidates, 
are responsible for the origin of the excess positron and electron fluxes and the observation of the 
DAMA annually modulated signal. It is also shown that the DAMA results and the other null 
results from direct searches for dark matter can be reconciled if 3 GeV < m^> < 8 GeV. On the 
other hand, in addition to SU(2)l doublet Higgs field H, the light singlet scalar field $, which 
is demanded to successfully construct the coexisting two-particle dark matter scenario and whose 
mass is taken to be just below 1 GeV, may play an essential role in resolving the PAMELA, ATIC 
and/or Fermi-LAT anomalies. 

PACS numbers: 98.80.-k, 95.35.+d, 14.60.St, 14.80.Cp 
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[l] Introduction. The quest for identification of the 
missing mass of our universe is one of the most fundamen- 
tal issue in astroparticle physics and cosmology. The evi- 
dence for non-baryonic dark matter (DM) inferred from a 
combination of cosmological and astrophysical phenom- 
ena becomes more and more convincing, which alludes 
the existence of new physics beyond the standard model 
(SM) . Recently, the DAMA collaboration reported strong 
evidence for model independent annual modulation sig- 
nature in the rate of single scattering events by combin- 
ing data from DAMA/Nal and DAMA/Libra with 8.2a 
significance [l[. This modulation signal can be inter- 
preted by the scattering of weakly interacting massive 
particles (WIMPs) forming the DM halo of our galaxy. 
While the interpretation of the DAMA results in terms 
of the DM-nuclcon elastic scattering is almost in con- 
flict with the null results of various other direct detec- 
tion experiments 0, Q , there is a small region of WIMP 
parameter space which can simultaneously accommodate 
DAMA experiment and the null results of the other di- 
rect DM detection experiments [J, Q . On the other hand, 
recent experiments measuring high energy cosmic rays 
may offer important insights into the existence of DM. 
The PAMELA experiment has presented a significant 
positron excess flux over the expected background with 
no excess in the corresponding anti-proton flux [6j. This 
result confirmed previous observation of a positron excess 
by HEAT and AMS 0. The ATIC/PPB-BETS experi- 
ment has shown significant excess electron and positron 
flux at energies around 300-800 GeV [1, ©]. More re- 
cently, Fermi-LAT experiment have also shown an excess 
electron and positron flux in the same energy range as 
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in ATIC but its strength was not strong compared to 
ATIC [13 ■ Since the positron excess observed in the 
PAMELA experiment is not accompanied by an anti- 
proton excess, an unified understanding of those exper- 
imental results based on DM demands new interactions 
through which the DM primarily self annihilates [n], EH 
(or decay [H, 0) into only leptons, and a 'boost' factor 
which could either have an astrophysical origin or have 
a particle physics origin such as Sommerfeld enhance- 
ment because the required annihilation cross section in 
the galactic halo is orders of magnitude larger than the 
thermal relic expectation [Til. [l5lTl^|. So, it is likely that 
the DAMA annually modulated signal is not reconciled 
with the cosmic ray positron and electron excess in the 
framework of one and only one DM scenarios. 

Recently, we have proposed an extended seesaw model 
to simultaneously and naturally accommodate tiny neu- 
trino masses, low scale leptogenesis and dark matter can- 
didate by introducing extra singlet neutrinos and singlet 
scal ar p articles on top of the canonical seesaw model 
[T7I filf . Furthermore, we have proposed a doubly co- 
existing two-particle DM scenario [lj| by allowing both 
an extra singlet Majorana neutrino S and a light sin- 
glet scalar particle ip two DM candidates. Such a sce- 
nario containing more than one DM may be desirable in 
the case that there exist a few incompatible phenomena 
which are very hard to reconcile in the scenarios with 
only one DM. 

The purpose of this letter is to investigate how both 
DAMA experimental results and the cosmic ray positron 
and electron excess observed from PAMELA, ATIC 
and/or Fermi-LAT experiments are simultaneously ex- 
plained in the extended seesaw model with doubly coex- 
isting two-particle DM proposed in [l{J. In this work, we 
slightly modify the model proposed in [l^ | by replacing 
extra singlet Majorana neutrino with singlet vector-like 
neutrinos so as to simply resolve the cosmic ray anomaly 
while keeping to accommodate tiny neutrino masses and 
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low scale leptogenesis. We notice that to achieve our 
doubly coexisting two-particle DM scenario in the renor- 
malizable framework as shown in [l!|, an extra singlet 
Higgs scalar field $ is necessarily introduced, which may 
open up new channels of DM annihilations. As will be 
shown later, in this scenario, this scalar field $ may 
play an essential role in resolving the unexpected electron 
and positron fluxes measured at PAMELA, ATIC and/or 
Fermi-LAT if the mass of <& has rather small around just 
below 1 GeV so as for the annihilation cross section to be 
enhanced via a mechanism first described by Sommerfeld 
0) EE EH- Once this new force carrier $ is included, 
the possibility of a new dominant annihilation channel 
SS — > opens up. The $ mixes with the Higgs allow- 
ing it to decay into the final state fermions, and if the <!> is 
taken to be light, it is kinematically constrained to decay 
to mostly lepton pairs preventing from producing anti- 
protons, so that the excess of positron and/or electron 
observed can be accounted for. In addition, the DAMA 
results will be explained by considering the singlet scalar 
t[> as the lightest DM candidate with mass of order a few 
GeV. Thus, the DAMA results, the excess positron and 
electron fluxes produced from the cosmic rays, low scale 
leptogenesis and light neutrino masses can be simultane- 
ously accommodated in our model proposed. 

The doubly coexisting two-particle DM scenario is 
achieved by the following Lagrangian 



C 



C + (Y D LHN + Y s NiPS + h.c.) + M N N T N 



1 



X, 



Y^SQS -moSS+ -ml ip 2 - -^ip 4 - XH^Hip 



~m| $ 2 -^$ 4 -A 3 V 2 $ 2 



A_2. 

4 



XaH^H® 2 



(1) 



where the first term is the Lagrangian of the SM and ki- 
netic terms of the singlet fields, and L, N, S and ip stand 
for SU(2)l lepton doublet, singlet heavy Majorana neu- 
trino, singlet vector-like neutrino composed of two Weyl 
fermions, and light singlet scalar, respectively. Note that 
S and ip are our doubly coexisting dark matter candi- 
dates. Finally H and <& denote the SU(2)l doublet and 
singlet (Higgs) scalar fields, and m$ is assumed to be 
smaller than 1 GeV. 

In this model, light neutrino masses can be achieved by 
the usual canonical seesaw mechanism and low scale lep- 
togenesis of order 1-10 TeV leptogenesis can be realized 
via the decay of the lightest heav y M ajorana neutrino N. 
Similar to the case in Ref. (l7l Il8j. there exists a new 
contribution to lepton asymmetry e\ arisen from a self en- 
ergy correction of the vertex generated due to the Yukawa 
interaction term Y$NipS. For Mn x ~ M^ 2 < Mjv 3 , the 
lepton asymmetry is approximately given by 
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Since the size of (Ys) 2 i is not constrained by the out-of- 
equilibrium condition, large value of (Yg) 2i is allowed for 
which the second term of Eq. (|2|) can dominate over the 
first one and thus the size of E\ can be enhanced. Again 
similar to Ref. [13, EH , the successful leptogenesis can 
be achieved for Mn 1 of order a few TeV, provided that 

(KsW(vs)m ~ 103 and M k/ M k ~ io. 

In order to guarantee the stability of the 2DM candi- 
dates, we impose the discrete symmetry Z 2 x Z' 2 under 
which all the SM particles and $ are (+, +), the singlet 
neutrino S is (— ,+) and the singlet scalar boson ip is 
(+, — ). Now, we demand that the minimum of the scalar 
potential is bounded from below so as to guarantee the 
existence of vacuum and the minimum of the scalar po- 
tential must spontaneously break the electroweak gauge 
group, <i/°>,<$>^0, but must not break Z 2 x Z 2 
symmetry imposed above. 

After spontaneous symmetry breaking, the part of the 
scalar potential is given by 

V = ^m 2 b 'ip 2 - \^n\h 2 - ^m 2 2 + 2X 4 v h v 4 ,h(p 
Xs ,a Xi 3 Ai 4 X 2 4 

+ + ~ Vhh +J6 h + T 

+ X 2 v <t> (p 3 + ^ 2 h 2 + Xv h hip 2 + X 3 ip 2 <p 2 + 2X 3 v tj) (pip 2 
+ ^-h 2 (p 2 + XiV^cp + X A v h h(p 2 + h.c, (3) 
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3A2^ — X4V 2 . Here, we have adopted 

\/2H T = (h, 0) and shifted the Higgs boson H and the 
singlet Higgs scalar $ by H — > h + Vh and + 
respectively. Since there exists a mixing mass term 
between h and (p, we rotate them with <p — sh' + ccp' 
and h = ch' — s(p', where s and c are sin# and cos 9, 
respectively. 

For rritf, < 1 GeV and m s >> m^, the singlet neutrinos 
S annihilate into mostly <p<p due to the cubic coupling of (p 
in Eq. ((3]). The </>'s can then subsequently decay into SM 
particles, which arises due to their mixing with the Higgs 
field h. For the case of = 0.25 GeV, the (p mostly de- 
cays to muon pairs, which in turn produce electrons and 
positrons, and thus the resulting spectra for the electrons 
and positrons are much harder than typical e + e~ spectra 
co ming from weak-scale WIMP annihilation as shown in 

HEEl. 

The amount of cold dark matter in the Universe, which 
has been determined precisely from 5 year WMAP data 
[H, is given by flcDAih 2 = 0.1099 ± 0.0062. Assuming 
the coexistence of two dark matter candidates, the relic 
abundance observed must be composed of the contribu- 
tions of both S and ip, ^ s ^ 2 + Sl^h 2 = f^cDM^ 2 - The 
relic density of each dark matter species is approximately 
given by (0.1p6)/ < av >j (i = s,ip), where 

< (tv >j is the thermally averaged product of its annihi- 
lation cross section with its velocity. For our convenience, 
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we define the parameter Si as a ratio of flih 2 to SIcdm^ 2 , 



Clih 2 



(4) 



where e s + = 1. In fact, the parameter Si represents 
the fraction of the mass density of each dark matter 
species in our local dark-matter halo as well as in 
the Universe. Since the values of Si are unknown, we 
consider a few cases by choosing their values in the 
analysis. Each f^ft 2 can be calculate with the help of 
the micrOMEGAs 2.0.7 program [2l| by taking input 
parameters appropriately. 

[2] Implication for DAM A experiment: In order to 
interpret DAMA results [l[ in terms of DM-nucleon scat- 
tering, we choose ip, the lighter DM particle of 2DM, 
to be relevant for the experiment. Note that the heav- 
ier DM S of order a few TeV is also demanded in order 
to explain the high energy cosmic ray anomalies later. 
To investigate the implication for DAMA experimental 
results, we first have to estimate the DM-nucleon elas- 
tic scattering cross section predicted in our scenario. So 
far most experimental limits of the direct detection have 
been given in terms of the scattering cross section per 
nucleon under the assumption that there exists only one 
DM candidate. In the scenario of 2DM, the cross section 
for the WIMP-nucleon elastic scattering a e i is composed 
of <r s and o~ [22l |; 
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with too being the WIMP mass, where we set = 
as the relevant DM mass for the DAMA experiment. 

In our model, the non-relativistic 5-nucleon elastic 
scattering cross section is given by 



<t s (nucleon) 



1 

-in 



sin 29Y t pm s m'^ l f 



(m„ + m s )vh 



-i 2 


" 1 


1 I 




4 






lm h 





> (6) 



where m n is a nucleon mass and / is defined by the rela- 
tion fm n =< n\ J2 q m g<rel n > whose size is determined 
by Hl|, 0.13 < / < 0.62. For our numerical estimate, 
we take / to be 0.36. The first and second terms in the 
parenthesis correspond to the elastic scattering mediated 
by the Higgs field h and scalar field cj>, respectively. In 
the case of scalar ^-nucleon elastic scattering, the non- 
relativistic elastic scattering cross section for ip is given 
by 
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where A' = XvhC + 2\^v^s and A" = — XvhS + 2X^06. 

In Fig. H]-(a), the pink-colored rectangular area 
presents the predicted region of the parameter space 
(a e i — m^) in our model for several fixed input pa- 
rameters given in the panel. Here, we restricted the 
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FIG. 1: (a)DM-nucleon elastic scattering cross section vs. DM 
mass. The DAMA experimental results are presented by the 
grey-colored regions and the pink-colored region corresponds 
to the prediction of our scenario for given input values pre- 
sented in the panel and 3 GeV < < 8 GeV. (b) Allowed 
regions of the parameter space (tan 8, A) from the fit to the 
DAMA results combined with the other null experiments for 
=3, 6, 8 GeV and the same input parameters as in (a). 



region of to be 3 GeV < m.0 < 8 GeV. We see 
that there exist some parameter regions compatible 
with DAMA experimental results as well as other null 
experimental results. Fig. [l]-(b) represent the allowed 
regions of the parameter space (tan 9, A) from the fit 
to the DAMA results combined with the other null 
direct search experiments for m^p =3, 6, 8 GeV and 
the same input parameters as in Fig. [T]-(a). Before 
discussing about the indirect experiments on DM, 
it is worthwhile to notice that such a small mass 
of ip taken to be consistent with the DAMA results 
does not affect e~ /e + fluxes observed in the galactic halo. 

[3]Positron excess from DM annihilations: Now, let 
us show that the PAMELA, ATIC and Fermi-LAT data 
can be accounted for by regarding singlet fermion S as 
a relevant dark matter much heavier than ip, which an- 
nihilates into dominantly cfxp, and then the 0's subse- 
quently decay into mostly [i^ [i~ when is taken to 
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be 0.25 GeV. In order to calculate the galactic cosmic 
ray (CR) propagation, we use GALPROP program [24| 
which simulates the propagation of both cosmic rays and 
DM annihilation products in the galaxy. The propaga- 
tion equation for all CR species is given in [24( . To solve 
the propagation equation under the assumption of free 
escape of particles at the halo boundaries, we fix the val- 
ues of the parameters in the propagation equation as fol- 
lows: the diffusion coefficient Dq = 6.10 x 10 28 cm 2 s _1 , 
6 = 0.33, the Alfven speed is 30 km s _1 , the convec- 
tion velocity is 5 km s _1 , the gradient of convection 
velocity is 7 km s^kpc™ 1 and the normalized electron 
flux is 1.325 x 10~ 8 cm" 2 sr- 1 s- 1 MeV" 1 at 10.6 GeV. 
We solve the propagation equation by using the code 
gladef_50p_599278 placed in GALPROP webpage [2j|. 
In addition we use an NFW density profile [26|], so that 
the core radius and the local DM density are taken to be 
20.0 kpc and 0.3 GeVcm -3 , respectively. 
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FIG. 2: (a) The ratio of positron to electron plus positron 
fluxes and (b) total electron plus positron fluxes, arising from 
the annihilations SS — > (f><p and then <f> — » fi + fi~ . B stands 
for the boost factor relative to < av >= 4.56 x 10 -26 cm 3 s _1 
which is satisfied thermal relic abundance for e s =0.9 



In Fig. [2j we present the predictions of our scenario for 
(a) the ratio of positron to electron plus positron fluxes 
and (b) the total electron plus positron fluxes, which are 
originated from SS — > <jxj>, and subsequent decays of the 
f/Vs into fi [i~ for the same input values of the param- 
eters. As for input values, we take e s = 0.9, m s = 2.5 
TeV, = 0.25 GeV, and < av >= 4.56 x lCT^cro^ -1 
which satisfies the thermal relic density of S for the given 
e s . Then, the contribution of S to the local DM den- 
sity (p° ) is 0.27 GeV cm" 3 while that of ip to the local 
DM density is 0.03 GeV cirr 3 . In these estimates, 
we invoke the boost factor (£>) reflecting Sommerfeld en- 
hancement through which the halo annihilation rate is 
enhanced, and it is given by 
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where a lies between 10 3 and 10 1 [16[. The each 
curve in Fig. [2] corresponds to different boost factor, 
B=1535 (1154, 640) for black solid (red dashed, blue 
solid) curve. The red dots with error bar correspond 
to the measurements from the PAMELA (Fig. [2r( a ))- 
The pink (grey), blue (black) circles and black crosses 
in Fig. H]-(b) correspond to the measurements from 
ATIC, PPB-BETS and Fermi-LAT, respectively. As 
one can see from Fig. [U-(a), the PAMELA data for 
the positron fraction is consistent with the black solid 
curve corresponding to B=1535. For the same value B 
(=1535) as in Fig. d-(a), the prediction of E 3 dN/dE 
as a function of E appears to be consistent with the 
ATIC data as well as the PPB-BETS data, whereas it 
is inconsistent with the Fermi-LAT data. On the other 
hand, we see that the prediction with £>=640 (blue 
curve) is consistent with both ATIC data (except for the 
peak) and Fermi-LAT data, whereas it is inconsistent 
with the PAMELA data. Therefore, it looks rather 
difficult to accommodate the PAMELA, ATIC and 
Fermi-LAT data simultaneously, which may imply that 
there exist other astronomical sources [13, [28], [2t| |3(| • 

[A] In conclusion, we have shown that the extended 
seesaw model proposed to achieve low scale leptogenesis 
can resolve the anomalies in the indirect detections of 
annihilation products observed from PAMELA, ATIC 
and/or Fermi-LAT and simultaneously accommodate 
the recent DAM A experimental results. In this model, 
an extra vector-like singlet neutrino S and a singlet 
light scalar field ifj which are coexisting dark matter 
candidates, are responsible for the origin of the excess 
positron and electron fluxes and the observation of the 
annual modulation signature at DAMA. Furthermore, it 
has been shown that the DAMA results and the other 
null results from direct searches for dark matter can be 
reconciled if 3 GeV < < 8 GeV. On the other hand, 
the light singlet Higgs scalar field <f>, which is demanded 
to successfully construct the coexisting 2-DM scenario 
and whose mass is taken to be less than 1 GeV, may 
play an essential role in resolving the PAMELA, ATIC 
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and/or Fermi-LAT anomalies. 
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